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ABSTRACT 


A software package for general analysis of real time 
series for geophysical data was developed. The package 
consists of FORTRAN IV callable subprograms that are as- 
meoned diiferent tasks such as input, output (either 
peamced Or graphical), initial data analysis and filter- 
mre, spurious data reyection, FFT, and spectral analysis. 
MeecOopnysical problem involving detection of unusual deep 


ocean pressures was analyzed. 
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Pane ODUCREON 


Lic ie Ome cimes series Gonsidered here are geo- 
mipsi cal pvocesses Usually, time 1s the independent vari- 
amen and such quantities as velocity and direction of flow, 
mempetature, and pressure are the dependent variables. In 
miles analysis and design of geophysical experiments, it is 
useful to have at one's disposal a set of software sub- 
momeames tO aidedn this analysis. Thus, a software pack- 
age for general analysis of real time series geophysical 
Mea was developed. tJhe package consists of FORTRAN IV 
@am@etable subprograms Dit@raremaisoiened ditterTentgtasks pouen 
Bemenipue, Oucput (Either printed @7 graphical) geinitial 
DE MEGiG SS anata -¢Crigic » SpuUnblous, cata @e,Cotlonm are es 
eeectral analysis, and so forth. A geophysical problem 
mime lving detection.of unusual deep O€ean Pressures was 


ame yzed. 





[i ner in 1 ol Cee RODEN 


A. Sie EMENT OhewnE ByROBLEM 

Ticoucheut mature Ehere exist many phenomenon that 
occur only as "unusual events." In the realm of oceano- 
REapny one such unusual event is called a tsunami. The 
word tsunami comes from the Japanesse and means storm wave. 
mEeis known that the presence of a tsunami is indicated by 
a definite change in water pressure. Other information 
momeectning the velocity and direction of the current, water 
momperature, salinity, and so forth, would also be useful 
in the analysis and understanding of these events. However, 
Mieeintrequent rate at which tsunamis occur makes their 
Bees vation difficult. In order to anel 
these events, methods were devised to detect and record 
miei Lhe problem then is to logically design a detector 
that will automatically sense the presence of a tsunami 


miomeause all the pertinent data to be recorded. 


B A MODEL 
mmonder tO determine the prasence -OLwaneunicuaimey Chi, 
a detector which receives digitized pressure data was de- 
waned sit |(e Throughout the following discussion rerer to 
Biwi Mereneiic Mocmkemaecioneet =the detector... |herde— 
Pe COLucanmoe Civ Gedmineousuxalogteally Separate Sections. 
ate titmniehe em resSULC =sensSOns OL Vani me 


deuce co somes OMn St leatnhom ond reliapr lity are presently 
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available on the commercial market. Digital sensors 
ma-meralivecount the number of cycles of a pressure sensi- 
tive oscillator over a given time interval. This summa- 
momen C€frect integrates the pressure signal. fhe 
mmpecration process serves as a low pass filter. That 
is, high frequency pressure variations relative to the 
Eaupling interval tend to average out, thus contributing 
equally to each sample. This filtering of high frequency 
feiartons Might prove undesirable in some applications. 
es this particular problem, however, almost all of the 
ieen L£requency variations are naturally filtered out be- 
@ause Of Sr enidtion Of WReSSUmenVarlaulons dlc .O sO cea 
maepth. 

it was determined that the frequency band in wnich 
Beuinmanis Occur extends from % ne SUNGyeGles sped wma, 
Nyquist frequency theory indicates that in order to detect 
ferecauency Of N cph, a sampling rate of at least ZN cph 
is required. Hence the sampling interval was determined 
emee One Sample per minute. 

PiplweasuEber. “An impwt bUftcmeneceives the most 1e- 
Cent pressure reading in digitized form. The actual hard- 
ware implementation of the input buffer 1S oO Cee kone eine 
problem. Conceivably the input buffer could be anything 
moom am acoustic delay line to a magnetic core. 

The buffer is large enough to retain the most recent 
121 consecutive samples. Hence, a two hour history of 


pressure data is available for analysis at any one tine. 





Peoure 2(a) 1S a plot of typical pressure data over an 80 
molt period. 

Pe oee eer Mele mm imcmprcociee Of ay tsunami could 
be indicated by an abnormal change in pressure from one 
Boame in time to another. The phrase "change in pressure" 
indicates that the derivative should be computed. As a 
Mose eapproximation to the derivative, the positive first 
difference was determined. 

MmacorpoOraced wichim the first differencer was a spur- 
mous data gate. Here each first difference was compared 
mast a preset threshold value. If the threshold value 
was exceeded, then a sensor error (spurious error) was the 
muopaple Cause since deep-ocean pressure changes are rela- 


+aaral 
eww we sh 


aoe Lon lt ees oe Te SUC ait Cimomcens 2125 e 
difference might falsely indicate the presence of a tsunami, 
am excessive first difference was set equal to the previous 
mmeseeaiiterence. The threshold value was determined after 
analysis of test data. First differencing, as computing 
Mme tirst derivative, removes the mean pressure value and 
mmesses drift. Fissure Z(b) illustrates the effect of 
meece differencing on the data of Figure 2(a). 

ace ites. eee DOS cmOnmmEnic mt Fle tite i ts se 
memovel tne Effects of the diurnal and semi-diurnal tides. 
This was accomplished by forming a weighted sum from the 
Sample data at 4 hour intervals. The coefficients were 


determined from analysis of the tidal spectrum. The tide 


miveer OUuLpUG EMen 1S : 
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D, = dP, - 3.920392+dP,_,, + 5.841800-dP,_ 


30 60 


zs 9 2 Dp 
920597 "dP on, * 8P i 45) 


where dP. ime brecteGiikerence at tame t {2 | . 

Line cC)mgnewoGicecatectmor the tide filter on 
miteetirst differenced data of Figure 2{b). Note how the 
tide filter extended the effects of a spurious data value 
Me there was no spurious data gating. 

Reetiracatton., As Preure 2Z(c} shows, Ehe output of 
@mes tide filter presents an erratic picture of the changes 
in pressure. In order to observe any trends in this Sco 
an integral over some fixed period of time needed to be 
@emputed. Any integration of the tide filtered data 
mela only produce a result close to zero. This problem 
was overcome by taking the absolute value of the data in 
two ways; that is, the negative data values were made 
positive so that the data was no longer distributed about 
zero. The most straightforward method was to simply ignore 
mecren (1.6€., rectification) of each number outputted 
meen the Made filter. The second method was to square each 
weve. Figure 2(d) shows che Quit Of SGUariIne= the frde 
mere r OULDUT - 

Integration and Detection. The rectified or squared 
tide filtered data was then integrated. The method of 
integration was based on the solution of a differential 


uct ion slice eencrial tTormyor the Integrating function 1s: 


ml 





Vee so enty Oe) wt tamer it, 25) ceo MINUS 


I 
= 


and where yy 


a= a time constant associated with the inter- 
Vici cam Ge ood eaceme, 
x, = the data to: be integrated (rectified or 
Squianed OUEDUtGmes seEnae 1 lter)= 
The value of a was determined experimentally via the simu- 
lation (a = 1/60). Figure 2(e) is a plot of the integrated 
imme tiitered pressure data after the squaring process. 
Dependent On the value Of of, od Grmeereal threshold 
Mmememiead co be chosen such that whenever the integrated 
fmeenal exceeded this value, a switch was closed causing 
"the contents of the buffer to be recorded on a tape or 
Qumerwise be permaneintiy recorded. At the end Gt a two 
hour period of normal input the buffer would once again 
Mmememccordacd. At this point the value of the integrated 
BSmomal would be tested again. If the integrated signal 
Prmeeectil!l greater than the threshold value, then the butr- 
mer wOUld again be recorded after two hours; otherwise 
meter cwo hours the testing is resumed as before. This 
Beechmigue guarantees that a minimum four hour record will 


be available for each event detected. 
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PeWhEe2. (alePressure, (b) First Ditferenced Pressure (spur; 


TousmiebeaGate. Off) s(c) Tide faltered Pressure, 
(dyyiSawered Tide Filtered Data, (e} Integrated 
Slcidemyomsus Time sin Hours. 
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III. A SOFTWARE SIMUCATION 


ie GENERAL APPROACH 

To be generally useful a software package must be able 
Memcimulate a real system on the level of logical components. 
ite 1S, the basic information extracting procedures of the 
system must be available for graphical or statistical anal- 
ysis during the simulation. Accordingly, the following set 
@mecaehteen FORTRAN callable subprograms was developed. 
The subprograms and their associated formal parameter are 
listed alphabetically in Table I. The subprogram UNIFRM 
2S a function subprogram and XYAXIS is an entry in DRAW1. 
Sieterwise all the subprograms listed in Table 1 are sub- 
Meutines, The various subprograms are explained in the 
foerowing Sections according to their use in the unusual 
wees problem. Each subprogram is listed in its entirety 
mieeeme Appendix. 

bach SuUbrOWwEame wses the Varreble FILE as ats first 
formal parameter. The meaning and use of FILE will be 
Pues Once here to avoid duplication in the Suck Tome which 
follow. 

iets wierd lovaived VYartaple tdentifying tie secur Teme 
wea tile. This number is for user identification purposes 
eee tubes usually set equal to 1.0 or the DO-loop in- 
@ex if the simulation is being performed iteratively. When 
a subyoutine iS called, a message is printed indicating 
that the Subroutine has completed its task. The message is 


Oueeene 16 EM 
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SUBROUTINE 


ANTIDE 
CARDS 
DRAWIL 
me REAL 
oem 
Eee | D 
HARMON 
INTGRL 
surPurT 


PACK 


POWER 
RECSQR 
RESET 
SCAN 


StAlS 


TAPE 
UNI FRM 


XYAXIS 





FORMAL PARAMETERS PAGE 
CRI, IU Cpa 1h) 26 
(FILE, Y, L) 16 





CUE ee Ne LOG pie X11)! 122 












































(FILE, A, B, N, ISN) 23 
(FILE, A, B, N, ISN) 23 
(FILE, Y, DIFF, L, THRESH) 25 
CETL. Yo ioe Leelee Hey) 24 
GBR CC jbeeieuete D7 
(FILE, Y, L) 28 
(EDbES) Yeey 2) 2 DEN se RAMP 

ISEED) Lis 
GHUWE, SO ins Cas) cl) 25 
OUD vs Ve ig SY 26 
(GE TIDE ee Yaxcila aN 2 eles) 18 
(FILE, Y, L, YMEAN, YSTDEV, L2) _ Li? 
(FILE, Y, L, YMAX, YMIN, YMEAN, 

YSTDEV) ey 
GRIEG VS IL mG 
(ISEED) : ay, 
CONES SONG | 29 





WADIEE VAL: 


Subroutine Names, Formal Parameters, and Page References. 


iS 





PubrEe NOeeXA= ks message’. 


Moe digits XX.X indicate the current value of FILE. Execu- 


mon Of SUbroutine DRAW! causes FILE to be incremented by 


i. 1. 


BD. DATA INPUT 


CALL CARDS (FILE, Y, L) 


CALL TAPE 


(FILE, Y, L) 


Data inpuc 1S accomplished by ecallinovene Orethe above 


subroutines. 
s.: 
Rael ec - 
Y 2 
L = 


The meaning of each of the formal parameters 


Soom ocm@t 1 Omen! | sae 

reat esting ly ad imens honed wag t ayy sis sonny ileln 
ine wildted ells. “wemibemnaciid: 

Ae Nie Cee Wahue devant pilkem oa Ud ln lemme le sone 
Sie SNORE, Leh subroutine CARBS, oo MUS bac auamn 
the NUMpemw Ot <datanec lemen £5.) Os Casteera ci meie mein 
SUbTOULIMe TAPE YL may sequal eeneenumD eines 
diate aie ee men: S) (tO ber C'al clas iehly mrOnigumlat ihc equal 
them mMaxd mlm) Sieze sO Y Sas edetaiimed iis c ic mei Ne 
STON statement of the main program. If an 
end-of-file mark is encountered on the tape 
before the L data elements have been read, the 
SUD GOMeINeG Nadees aAncerelUYTis. Eile. ume deme 6 cilC. 


NoMmecmictmolilvenrocessed = as the value of ile 


inesc subrouUGames shoulld be reloaded with the proper 


DORMbele Stabemcnts Inserted becrore the simulation 1s run. 


PEnenehy seme wrOkKMAT as (10F 7.041020) - 


16 





oF INITIAL SCAN OF DATA 


CALL STATS (FILE, Y, L, YMAX, YMIN, YMEAN, YSTDEV) 


CALL SCAN 


(FILE, Y, L, YMEAN, YSTDEV, L2) 


Pe eis SOMeLamMesmtne = ecdscuwimcnescOpIySieGal data as being 


meeorded that the initial start up of the sensing device 


intorudces spurious data which is recorded as valid data. 


Used in series, subroutines STATS and SCAN can determine 


Mrmextent Of any start Upp moise present in the data. 


iitcamcanine Of Cac of the tennal pa;anetersm is 


2 ie! 


Y 


YMAX 


YMIN 


YMEAN 


YSTDEV 


SIS eee areye Vt 

eed vaeltied wiSmnCiymcdinen SmOneG iammanmeeOils 
Edie tie datd sto Dem@dmal yzcd nice tiles 
subroutine alters the contents of array Y in 
any way. 

AM Antec er-Vallied sVararadbie equate O | time smEac 
of array Y. 

A real-valued nuinber returned by sub out ime 
STATS equal to the maximum value encountered 
I Bleie@ely 

A real-valued number returned by subroutine 
STATS equal to the minimum value encountered 
iiearray 

A real-valued number returned by subroutine 
STAUS equalevte themstatistical mean ot ene 
data in array Y. 

A real-valued number returned by subroutine 
SIMUS eaual to the standard “deviation Of the 


Gawain array. 
17 
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L2 - An integer-valued number returned from sub- 
PemrEinewocAN equal to the number of noise 
dependent elements in array Y as determined 
Peon tie fLollowine procedure. 

1) set I=l, set ITEST=0. 

2) ENE S re ieee SI DEY then set 
eroW 0 us crt. ela Tes 
otherwise, set ITEST=ITEST#+1. 

See LEM ESole -  LOCREneCne return. 6 elem naece 


Seu =] ecommerce) ae. 


1D. ADJUSTING THE SAMPLING INTERVAL 

CALL RESET Geo Seem Coasts arn 7 tu 

xvii PACK GIR A Oe iA emg lec Booed rege! 210) 

UNI FRM CL SISTED) 

Mine Semple size Of data EOmbewand) seamen pe weaagtis ted 
in two basic ways. First, the sampling interval may be 
Samieeed £LOrward by any specified amount. For example, if 
M@eewaminary analysis of the data indicates that the initial 
fea Clements are in error, then the data file can be re- 
fees) tliat the record begins latemuan the sequence there- 
by skipping the bad data. This 1s accomplished by a call 
Sresubroutine RESEI. ‘The LOVNalepaGame lel s (Of J bok lace 

IEE EIs Sess Seegeen, IML 

nf —Jieeedl- walued =) simea iy dimensmoncdear dy scon- 

Pine menenGata tO De reser. 
Ik =P ieee Oe GV allied variable Conall em ine si Ze 


Slwarndy.s \ . 
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189 ea St 


3s ; 


L2 : 


A real-valued, singly dimensioned array into 
which the reset elements of array Y are to be 
eaten AYrays Y amd Y2 may be the same 
physical array in which case the skipped data 
is overlayed and lost to further processing. 
An@anteger-valued variable. On the initial 
€qaee to RESET, LZ cquals the number ot data 
elements in array Y @ehich are to be skipped: 


Chol nec reba | 1b  erepbpellic. ilies ele elimechapen 1 2 


The second method of adjusting the sampling interval 


involves increasing the density of the data points. This 


1s accomplished by deal Om Subiaome Mice Gah Cm Onemam 


meamameters of 


2 = 


L2 : 


AO ho ager 
see Section TII.A. 

A real-valued. singly dimensioned anigay seo, 
taining the data to be packed. 

Ape einee cc fr - V all vie dea a teine ble me Cellet @ mm meee 
Onmar Cave oe 

A real-valued, singly dimensioned array into 
Wiehe Chie rexes andeds paleo ch ard acy Osea isiecy ann 
1s to be placed. ArraysyyY and 2 may pewrne 
same physical array. 

Ameimtecer valued Variable. -On. Giemwin rt ra 
COHEO Eile Sipreuttine, LZ must cata | aume 
maximum size of array Y2 as declared in the 
dimension statement of the main program. On 


recuuieeie will eotlalytre=size Of arrdy 12 


us 





as declared in the dimension statement of 
icmp EOchamem Om retumn, L2 will equal 
Chics Oren ay 9) cmt tne SUDTOULInNE Com 
Pee toCmolGceeasctiiive mw Ne sSmze O41 array YZ 
Calle Dene Corinimned from the sommula (L-1)x 
DOE Na ie 

IDEN Alive Vee) ule tums C Speci yimld etic 
NUE emo ec Mellen eNS hty timaturs eae 
Sired, The value ot IDEN in part determines 
Cie auicwon Ib Ze Cid toe CUmMmed tO. eme madam 
program. IDEN can also be interpreted to 
mean that (EN=-f points ane embedded petyeen 
Sach Ot thie rOrie tial polit o wan a1 ray 1. 

RAMP cee 2 lev ol ed ie et ab Ve sea vial hoy tien oa 


variability allowed when the embedded points 


0.0 then the 


B jee) lors alma one luce 5 Ilse jeunnyile 
embeddeds points »ame dete minpiecs oye ae cian 
he EDO de loOnme Del Weel suGGCS Sa \ Canny aumecs 
otherwise, the value ot RAMP is multiplied 
by a uniformly random number between -* and 
“my, and this value 15 added Co the interpolated 
Valliic 

ito I i ieinbecer - Vvaluled=Vvarranle titat Servesmas tne 
PANCdOMmMiUNDehesScea stom na random number Cenen 
Omer kl) SiO mie S x dae) ten begem 
equal to +3(mod 8). 

Conceptually subroutine PACK embeds the data as shown 

Bree doure 5. 
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pat. 


Y 
(PRESSURE) 


Y(I+1) Pi 





A 
Y (1) re 
IDEN-1 
embedded 
points 
I ea 
INDEX 


IG WUGRC el Shy 


Pemecptual Operation ot Subroutine PACK where RAMP > 0 
and IDEN = 4. 


ZA 





Function subprogram UNIFRM returns a uniformly random 
number between -¥% and 4%. The use cf UNIFRM is not re- 
Peametea tO Ssubroutime PACK, 

The original data used in the simulation was based on 
a four minute sampling interval. The detector, however, 
was designed to operate on one minute data. Hence, to ob- 
mom Ene proper Ssamplang interval the following call was 


mede on subroutine PACK. 
CONT IS RANCH Cy 2S IGG) er) (es al it ees re US, | CL Ot Jemma (0) oye elm 


LMAX equals the maximum size of array Y as dimensioned in 


Mmiem malin program. Note RAMP equals zero. 


-E. FOURIER ANALYSIS 
CALL FFREAL (FILE, A, B, N, ISN) 
CALL FFT (FILE, A, B, N, ISN) 


CALL HARMON (a OO Series ae he Go Paes eae el gee 2) 

ene POWER CED: Ys gi 

Three areas of Fourier analysis are provided for in 
the software package. Subroutine FFREAL returns the Fourier 
cosine and sine coefficients SETtte data Passed st Orie. oub. 
routine HARMON returns the approximation to the data passed 
moemmteon the basis of the Hl through H2 harmonics. Swub- 
memrine POWER returns the power spectrum Of the data. Sub- 
routine HARMON and POWER both call on FFREAL, and FFREAL in 
mr Galls On subroutine FRI. Subroutines FFREAL and FFT 
are variations of ea ae PSL Alva ore eis stee Weel Pieeyo mt 9 


mace router transcorm. 


Ze 





ine Lormals parameters of subroutine FFI are: 


Jee Ee : 


ia 8. = 


ISN ; 


Bee Section II1.A. 

Reaul-yvawwed, samraly idimensioned arrays. Ar- 
rays A and B originally hold the real and 
imaginary parts of the data and return the 
real and imaginary Fourier coefficients. 
Mimeneege ta ecm eimiadolc COUdIaEtOmENc alum 
ber of real (or imaginary) data elements. 

Ine Ene mene t | Oia enilsmina jer an eN injec 

i riheewe waviailicCm ieee rm | en iiaen C00 Cc monr 
ISN determined the step size of data in ar- 
Rayo end) ems Ome a Shoal iede feminine s 

VIOUS CNUs eh celavey “(Cleve dke rs “ere shivers Ligeti sll I 


be peri cummed . 


iiiestOrnal parameters for subpemtine FPRREAL are: 


fell bys : 


See SeCtilone lla 


aoe Shea aed. SiO yese Imes 1emcc waited \7s cement dle 


Original data values are alternately stored 
in arrays A and B. On return, arrays A and 
by Joe Gl Tele EOS TNe Bunch SIS MOTE s ec 


Ele Cieih ella penne mie(he satel cas te@rc ii: 


N Se ieee Tale dmv amadeaDle eG Uae: m Ct MicmmmlniiS 


ISN = 


Dex wGimewnemenes aim eanGay Aan In eiheweonmcc 
of this paper, N=L/2. 

An integer-valued variable. The magnitude of 
ISNeGetermines the Step size within arravs A 


deer seeon hac 0 kee a0 talon is eee l 


~ £ ato 





2 eee 


each consecutive element of the array is used 
If Siemeonputatlons ir the Magnitude 1s 2, 
CVG wOcic Wace Mnenem Ss UScd@mmHence, 1 the 
dgtamis Golmecinedswaenan @ Ssimgle array wasay 
X, rather than two separate ones, the call 


to FFREAL would be 
CAI RE Rae Abana ermeee(27)) 4.” Nee 621m 


Pieces Hon Ot Glo \iesmile fata Ve.) Clem te Mest) 


verse transform is computed. 


The formal parameters of subroutine HARMON are: 


BILE 


ye 


YH 


H1, H2 


SCCZOCC EO mums. 


A real-valued, singly dimensioned array orig- 


a+ oe ee ee 4-L leer 
AL L 


dere nas) (ae | ~~ mm ae vT FT 
Liane y CUILCaALMiite wie ee ariel is ce 


U se UdlLaus y Oo Ce 
On VetUyme arta, a) "COM cusniSeticm oumalc beer 
Sine and sine coefficients as returned from 


Joule) eglaA Lae 


- A real-valued, singly dimensioned array which 


on veturn holds the a@ppreximativon baged on 
tie hel sand tic. a pMOMGdhe 1 O1 see mcha wim O, bcslbnicalem, 
CONnmpeamnedmn etGay ee Ths  enOeccos milion Wes 
aedoup les Summadtlon. and could pe. Cine seons tile 
ing depending on the values of Hl and HZ. 

Ail ee en aeeuule de vaga ain C (OGta labo usGiMc monezc 
Oise Elicia VOM ojss larvae Jaye 

Positive integer-valued variables indicating 
Diidheimie CMOnlCsmetre ren rough H2e ine luisiue noe 
Gehaye’ ares tombce nse aved. 
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The formal parameters of subroutine POWER are: 

fea lee swe carpe Gr lem ill ae 

a Pee Ge man cuca mc Mig hy ScinenisLoned array eon- 
taining the data which is to have its power 
Spc Cea acOlpucd . POM return ws tne spower 
Spee enlm@mcoe |i lemen@orma Ine 911 |) tinecemaanst N—ily 2 
elements of array Y. 

L =e ees Chanelle @ mizcuelealb ne CCU nom &Mems sac 


Civ arraye . 


Fe FIRST DIFFERENCING 

CAE FIRSTD @SIOEIE Gy me OND cules Mei, ML BOR eS) cll) 

Bomourcince FIKSIDPecomputesmenc pOStelve stirs tec it 
ference of the data in array Y, and returns the array of 
ieee ditferences in array DIFF. The first and last ele= 
Mem@eseoL DIFF are zero. Subroutine FIRSID 1s so coded that 
mieeays Y and DIFF may Dew ehewsane physiC alan ra\e ns 6 yeesic ee 
mmomuURESH to any computed value, any first difference 
Meeeecr than THRESH will cause that first ditference to be 
taeeecqual to the previous first difference. 

ie formal parameters Of subroutine FIRSI) are: 

lo JL ibye: = see Secu Lon ti eae 

Y Awe a -Vakuecd as oly Cine ms HOMe Gama kane On 

anode Ca tc EOmbiem@—iadt Studd i beneneca: 

DIFF PevVercdlenvalvucd ciclo eames ONnedsarray wile 

Romo Ne pOstiivestirse ditferences,or the 


data in array Y. 
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ay 





L - An integer-valued variable equal to the size 
Stuimnay 1 (or array DIFF). 

THRESH - A real-valued variable against which each 
Pita POLTenee TS aconpanred..elt the current 
first difference is greater than THRESH 
(threshold value), cond the first difference 


Sees ee Ud leo mae mn Ige VOlOlLS. telieS tancr ite te miger 


G. POE PILTER 

CAEL ANT IDE (Ce ene | Dineen) 

Geopnysicalrevents Usually occur im the domain) of other 
meeoeinceresting phenomenon. In the case of the unusual 
eet agetector, one such phenomenon is the effect of the 
tides, both diurnal and semi-diurnal. Subroutine ANTIDE 
(anti-tide) filters out the combined effect of these tides. 

The formal parameters of subroutine ANTIDE are: 

RE = See peerlon. bl tA. 

i Ae red Varied, Si nel aime tse ine a scined ene Oils 
(CALI Tele Gee ee) seas Toile ed. 

EDIE 2 Ae oelay = Vase U.-aes aime y, jinene toned axeee into 
Winwe he cic st ideu if lbte hed dtd ins — ede ec Crm les 
rays Y and TIDE may be the same physical ar- 
ray. 

iL = AIemeececr-Valued variable equal to the size 


Cun yee Or dice Di) 


ale Ree EN Ge iteeslGNAL 
GENIGIES ROS KGIS ONE: GChiiiger Wo ee cwumbee LRG) 
Two methods of rectifying the data are provided. The 
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mirst anvolves taking the absolute yalue of each data value; 
wmmle the second squares each data value. 
Mice Ohidhwamaitecer Ss OLestmGouuemmne RrEcoQkR ({recti ty- 
sollare) are: 
PILE = Scewoce tno! Take 
ye Sheree valled as Limee a dt inensmomecd array Con- 
taining the data to be rectified or squared. 
ye - A real-valued, singly dimensioned array into 
which the rectified or squared data is to be 
placed. Arrays Y and YZ may be the same 
physical array. 
L “VAN Vin tever-Valuedyartaple sequal. «lon tie =steme 


Of alray Y (or arranges 


bipies - An intefer-Valued variatle €quai Co 1 or 2. 
Tf IRS*=]=. 1) then the data will be recti ied 
by Computime =the abso lute. vd lemme amme cieimc luc 
ment. if IRS = 2, then the data will be 
reCctitleds by ,GCOMputet ne seinem oCniaimeriG eae 
Che mive them 
I. INTEGRATION 


CALL INTGRL Dee eee, Js Ie) 

MitcemaelOnmiss Der rOnmnecd sn amu tap ly tno pies UNG £ 
MPeeevueusS data elements by a fading time constant (FIC); 
that 1S some positive number less than unity. The form of 


MicwiInteprating eCaquatLon 1S: 


VC. = EOD eee eS (SCR IS IA een ap 


oa 








where FTC is the fading time constant, and I is 


Ommrene current variable. 


Mice Ornampabaneters Ofesubroutine ENIGRL are: 


je dee = 


- ; 


aeG = 


fee OUTPUT 


See Sec iem Mi ke 

A real-valued, singly dimensioned array 
iron Celcom deammire ube min CCS raitoar 

A redali-y alvedees tiie. dimensioned array 
Wiel els anteeeceicecl Gece Ts jollece, 
ANePNLeCCeT=V a1 UCC mmataaable  COUAl tom chic 
Ou array X (On (arrays: 

A real-valued variable. <A sample value 


Re Cemilal be a / 6. 


ek OUTPUT (PILE, Ys 1) 


CALL DRAW1 


CALL XYAXIS  (XAL, YAL) 


imo Means Of CUTPpUuL are provided. 


Memes OUL the first L elements of array Y. The page 


the in@ex 


Con; 


Ineo 


s1ze 


of 


e | 
(FILE.9YouL, YX vit. LOG seche= on 


Subroutine OUTPUT 


for- 


treeeot OUTPUT is 200 numbers with index numbers per page. 


MicetoOnmal parameters of Subroutine OUTPU! Vare: 


elie eS Clemo el Onlmeiel lear 

4 ne Le dl Veet eS ile lye cd tiles NOMeGdmar hay sc omlt 
Pi omune Caea  tOwwe  primted. 

L noc to t= valle d ian lap le scdual CoOueneG. Milne 


Demme nc ONES OL Raunayerirenatl are v bo me 


Priced lay COU astire wt Ul bSizeror Farrday 


rimwinG CiieGcascwanlagcc Guant1ty OL Paper is 


BIG Ie LI eseigle@ehieiea! 
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SUD rOULIIC wn Ll OlOtSuene tigse Lb elements of array 


Y. The formal parameters of subroutine DRAWI are: 


jes! = 


y = 


YMAX, - 


YMIN 


LOG = 


mie Xl = 


sieve Sree eg IS 

A real-valued, singly dimensioned array con- 
taining the data to be plotted. Subroutine 
DRAW1 does not alter the contents of array Y. 
An integer-valued variable equal to the num- 
Bet 761 “data soot tcmeomy C70 Loe ted: 

Real-valued variables equal to the maximum 
AiG ene Vash US mom Cm MlOiet cl meme Laat ee 
iM Pee icien tic: NO feisSmateeisied led. 

AUG) Sugeno eMlvloel Wreldakelo be olla ee 0) ie SIL - 
lt LOG = Ty then the data in varray Y 1s plo. 


VC — 
oe) Vv 


ae 


BO Onl. ee 5 GoeeCe Olver ensic aynew 
tier d dtd Suelo t wediet mec leu ihe medic 1) eee 
Real-valued, simely Vdamens )onedm an Gay cmcon. 
tainine literal constants. the winaxsimumn seme 
String allowable is 40 Characters et seein 
plot label which may be blank> and Xloas the 


ser ss laene mine at Fone ieane le 


Nicmrchounss oO pene  plouraxccmane si lterdlvzcads toe) 


mienes £Or the x-axis and 8 inches for the y-axis. However, 


these dimensions may be changed by calling XYAXIS. The for- 


mal parameters of XYAXIS are: 


XAL = 


A real-valued variable equal to the new length 


Opie wxX-ax1S In -inenes- 
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YAL ee te ae Can ome dihtlmeo une =mew lene elt 
Ci sthe .) Cis mitietnehcocum line Value of  YAE 
SHmovule Ie less tee 9S maelies 
Subroutine DRAW] (and XYAXIS) make use of the standard 
pemet package for the CALCOMP 765 plotter available on the 


IBM 360 at the Naval Postgraduate School. 
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LV ee CONCED SEONS 
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Iie (ela setectad. iene question arose as to what was the 
best value of the fading time constant used in the integra- 
tion scheme. That is, given the form of the integrating 
equation and that integration is naturally a low pass fil- 
Mer, what value of the fading time constant would produce 
the most favorable signal-to-noise ratio? Three values of 
ie time constant were chosen, and the integration was per- 
mrmed three times on the squared output of the tide filter 
mmm a different time constant in effect Seven einelas INE 
value of the constant used to produce Figure 2 was 1/60. 
Bimeesipnificance of this value is that because of the struc- 
mime Of the integrating equation, the contribution of a 
seeele data value after one hour's integration would have 
meemereduced to one third of its original effect. This 
feeee, Of the time constant was used. The other two values 
fueeiene time constant considered were 1/30 and 1/120. 
maple JI shows the results of two signal-to-noise ratio 
computations using the three constants. The two hour 
"threshold" level to noise ratio was computed by finding a 
level such that the signal exceeded that level for exactly 

ive nours. 

The integrations were performed a second time but the 
mectitied rather than the squared output of the tide filter 
was used. The ratios are trom? thesemintcoratlons are 


meso civen im lable Il. 
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From examination of Table II, it was determined hat 
eevalue ot 1/60 for the fading time constant gave the most 
desirable results when the integration was performed on the 
m@mared Output of the tide filter. From inspection of the 
miepuce Statistics of the first differences, the standard 
Meweracion ot the tirst differenced pressure data was ap- 
Meeximately eight. Hence, a three Sigma value of 25 was 
chosen as the threshold value of the spurious data gate 
teenie srirst differencing subroutine. Figure 4 shows the 
mmcet Of the spurious data gate on the detector output 
when this threshold limit was set at 25. Figure 4 also. 
memreSents i oheaeee time constant value of 1/60 appa oe 
beethe integration of the squared output of the tide fil- 
mee ccMpare Fiewres 2 amd 4.) 

The critical threshold value of the unusual event de- 
Memon was determined to be three times that of the back- 
ground signal of Figure 4(e). Hence, the three events 
detected in Figure 4(e) will cause the recorder to be ac- 
mevatced. 

ies power Spectra 4] Ge the Cider tel ecme aie ecdeme Gusmao 
Outputs are shown on Figures 5 and 6, respectively. ee 
feure 5 the null power readings coincide with even cph, 
and the peaks coincide with odd cph. (ees: Si esiltesd cigees 
poem the characteristics of the tide filter. The nodes of 
the power spectral plot of the detector output do not ex- 
aetly acree withethose of the tide filter because of the 


omrect OF the intepration process. 
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FIGURE 4. (a) Pressure, (b) First Differenced Pressure (Spurious 
Data Gate On), (c) Tide Filtered Pressure, (d) Squared 
iptGewPmitercd Pressure. wic) Intesrated Pressure Versue 


Time in Hours. 
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